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We present here a scenario, based on a low reheating temperature Tr <C 100 MeV at the end of 
(the last episode of) inflation, in which the coupling of sterile neutrinos to active neutrinos can be as 
large as experimental bounds permit (thus making this neutrino "visible" in future experiments) . In 
previous models this coupling was forced to be very small to prevent a cosmological overabundance 
of sterile neutrinos. Here the abundance depends on how low the reheating temperature is. For 
example, the sterile neutrino required by the LSND result does not have any cosmological problem 
within our scenario. 

PACS numbers: 14.60.St, 98.80.Cq 



In inflationary models, the beginning of the radiation 
dominated era of the universe results from the decay of 
coherent oscillations of a scalar field and the subsequent 
thermalization of the decay products into a thermal bath 
with the so called "reheating temperature" Tr. This tem- 
perature may have been as low as 0.7 MeV [l| (a very re- 
cent analysis strengthens this bound to ~ 4 MeV Q). It 
is well known that a low reheating temperature inhibits 
the production of particles which are non-relativistic or 
decoupled at T ^ Tr 3]. The final number density 
of active neutrinos starts departing from the standard 
number for Tr < 8 MeV but stays within 10% of it for 
Tr > 5 MeV. For Tr = 1 MeV the number of tau- and 
muon- neutrinos is about 2.7% of the standard number. 
This would have allowed one of the active neutrinos to 
be a warm dark matter (WDM) candidate Q. Experi- 
mental bounds force now all three mostly-active neutrino 
masses to be in the range of hot dark matter (HDM), but 
we can use the same idea on sterile neutrinos. 

Sterile neutrinos without extra-standard model in- 
teractions are produced in the early universe through 
their mixing with active neutrinos (5J- Dodelson and 
Widrow [6( (see also Ref . [7j ) provided the first analytical 
calculation of the production of sterile neutrinos in the 
early universe, under the assumption (which we maintain 
here) of a negligible primordial lepton asymmetry. Fig. 2 
of Ref. [U shows that mostly-sterile neutrinos produced 
in this manner have an acceptable abundance only if their 
mixing with active neutrinos is very small, for example 
sin 2 28 < 10~ 7 for masses m s <; 1 keV. In the presence of 
a large lepton asymmetry, sterile neutrinos are produced 
resonantly 0, with a non-thermal spectrum which fa- 
vors low energies ("cool" dark matter candidate). 

The main idea of this letter is that the primordial 
abundance of sterile neutrinos does not necessarily im- 
pose their mixing to active neutrinos to be as small as 
usually believed. We can, thus, consider sterile neutrinos 
of any mass and coupling, as long as other experimental 
and cosmological bounds are satisfied. These neutrinos 
could, therefore, be revealed in future experiments. Only 



for simplicity we do not deal here with neutrinos heavier 
than 1 MeV or with a large chemical potential, or with 
reheating temperatures lower than 5 MeV, but our idea 
clearly applies to all of these cases too [l0| . 

In Ref. it is shown that most of the sterile 
neutrinos are produced at a temperature T max ~ 
133 (ms/keV) 1 / 3 MeV. Thus, if Tr < T max the produc- 
tion of sterile neutrinos is suppressed. We follow the cal- 
culations of Ref. [6j , but consider that the production of 
sterile neutrinos, through the conversion of active neutri- 
nos produced in collisions, starts when the temperature 
of the universe is Tr < T max . 

In the calculation the active neutrinos are assumed to 
have the usual thermal equilibrium distribution /a — 
(cxpE/T+ l) -1 , thus, following Ref. Q, we restrict our- 
selves to reheating temperatures Tr > 5 MeV. We also 
restrict ourselves to the case of sterile neutrinos with 
mass m s < 1 MeV, so that we do not need to consider 
their decays into electron pairs. These sterile neutrinos 
are, therefore, relativistic at production. 

In the approximation of two-neutrino mixing, sin# 
is the amplitude of the heavy mass eigenstate in the 
composition of the active neutrino flavor eigenstate v a , 
a = e,fi, t. 

For some range of masses which depend on Tr, one can 
neglect all matter effects, so the oscillations are as in the 
vacuum. For example, for Tr = 5 MeV, the specific value 
of the reheating temperature we use in this letter, this 
happens for m s > 0.2 eV (0.1 eV) for v e v s [y^r *-* 
v s ). In this case, the v s distribution function turns out 
to be 

f s (E,T) ~ 3.2 d a (j^f) sin 2 29 f a (E,T) 

(1) 

where d a — 1.13 for v a = v e and d a = 0.79 for v a = v^ T 
|ll| . This distribution results in a number fraction of 
sterile over active neutrinos plus antineutrinos 

f = 2* ~ io d a sin 2 20 V . (2) 
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FIG. 1: Bounds and sensitivity regions for v & «-> ^ s oscilla- 
tions. See text. 
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FIG. 2: Same as Fig. 1 for ^ MlT *-* za. For i/ T <-» /A the 
darkest gray-blue excluded region does not apply. See text. 



Notice that the number density of sterile neutrinos de- 
pends only on the active-sterile mixing angle and the re- 
heating temperature. A low reheating temperature in- 
sures a small sterile number density, even for very large 
active-sterile mixing angles, as large as other experi- 
mental bounds permit. This makes sterile neutrinos in 
our scenario potentially detectable in future experiments. 
The i^-number density is independent of the mass of 
the sterile neutrinos, contrary to the result of Ref. 
Thus, the mass density of non-relativistic sterile neutri- 
nos, Q, s h? = (m s n Va / p c )h 2 depends linearly on the mass 
and on sin 2 28, 




The condition Q s h 2 < Sljjj/fc 2 = 0.135 0] rejects the 
triangular dark gray region of masses and mixings shown 
in Figs. 1 and 2. The values of masses and mixings for 
which sterile neutrinos constitute 10% of the dark matter 
are also shown with a dotted line. 

Figs. 1 and 2 show bounds for v a = v e and v a — v^ T , 
respectively, and for Tr = 5 MeV. They show that v s 
in our scenario are viable HDM candidates, while neu- 
trinos with m s > 1 keV are disfavoured, if not rejected, 
as WDM or CDM, by bounds coming from supernovae 
cooling and astrophysical bounds due to radiative decays, 
explained below. 

Through v a <-> v s oscillations, sterile neutrinos can be 
produced in supernova (SN) cores and escape, carrying 
away a large amount of the released energy. The ob- 



servations of v e ,v e from SN1987A constrain the energy 
loss in v s and yield a bound on the mixing angle. For 
m s 45 keV, v a <-> v s oscillations are matter suppressed 
and the forbidden range is 0] 

0.22 keV < |m s (sin 2 2(9) 1/4 | < 17 keV . (4) 

For m s 45 keV, the matter effects are negligible and 
neutrinos oscillate as in vacuum. In this case, the forbid- 
den range is 0| 

7 x 10" 10 < sin 2 28 < 2 x 10~ 2 (5) 

These bounds exclude the diagonally hatched region with 
thin lines in Figs. 1 and 2. The effective matter poten- 
tial in the SN core for v e — > v s conversions, V m , might be 
driven to its zero equilibrium fixed point, V m ~ 0, during 
the explosion 0, . In this case the v a — » v s conver- 
sion happens as in vacuum. However, using the SN pa- 
rameters of Ref. 01 , this does not change the bound in 
Eqs.|3|andEl The SN1987A bound on neutrino radiative 
decays [l5|, excludes the region above the line labeled 
SMM (Solar Maximum Mission satellite) in both figures. 

If the sterile neutrinos produced in non-resonant v e — > 
v s conversion, in fact, carry away a sizable fraction of 
the energy emitted in a SN explosion, asymmetric emis- 
sion of v s due to the presence of a strong magnetic field, 
could explain the very large velocities of pulsars 0] (see 
diagonally hatched region with thick lines in Fig. 1). 

Having restricted ourselves to m s < 1 MeV, the dom- 
inant decay mode of the mostly-sterile v is into three 
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neutrinos. Assuming neutrinos are Majorana particles, 
a lifetime r equal to the lifetime of the universe tjj — 
4.32 x 10 17 s is indicated in the figures with the full thick 
line (for Dirac neutrinos, the lifetime is larger by a fac- 
tor of ~ 2). Equivalent lines corresponding to shorter 
or longer lifetimes can be easily obtained knowing that 
the lifetime is proportional to (sin 2 29 ml)^ 1 . The de- 
cay mode into a neutrino and a photon happens with a 
branching ratio 0.8 x 10~ 2 (this decay is not GIM sup- 
pressed, contrary to the decay of an active neutrino into 
another active neutrino and a photon ^(|). The diffuse 
extragalactic background radiation (DEBRA) imposes a 
bound |17| that, in the relevant range of masses, can be 
well approximated by 

L t < {E/0.05 McVy 1 (cm 2 sr s) _1 (6) 

where J 7 is the differential photon flux. 

For decays with r > t rec (t rec is the time of recombina- 
tion in the early universe; the line r = t rec is also shown 
in the figures) the bound obtained for unclustered neu- 
trinos would reject the region above the dot-dashed line 
labeled DEBRA in the figures. In particular for t > tjj, 
for unclustered sterile neutrinos, the bound would be 

( J^-) < o io <rv° ( ^XV /2 I (7 ) 

However, this bound affects neutrinos with mass m s <; 
100 eV which are gravitationally clustered. Therefore the 
bound is not entirely correct. But the actual estimate of 
how much of the diffuse photon background would be due 
to neutrinos which decay after structures in the universe 
form, is missing in the literature. 

Abazajian et al. ^| proposed to observe clusters of 
galaxies with the Chandra and XMM-Newton observa- 
tories, in their high sensitivity range for X-ray photon 
detection of 1-10 keV. They proposed to reach a detec- 
tion energy flux of 10~ 13 erg/(cm 2 s) with the Chandra 
observatory, which would allow to observe a monochro- 
matic signal from the Virgo cluster, if 

(i^)" 2 (^) ,/3 , (8 , 

for m s = 2-20 keV (horizontally hatched region with 
thin lines in the figures). Here the density fraction of 
sterile neutrinos within clusters is assumed to coincide 
with the cosmological energy fraction (TIs/VLdm)- 

The lack of distortions in the CMB spectrum due 
to neutrino radiative decays, excludes all the vertically 
hatched region with thick lines in the figures [19j . 

Structure formation arguments impose sterile neu- 
trinos which constitute the whole of the dark matter 
(WDM) to have m s > 2.9 keV jl. Note we use 2.9 keV 
instead of the 2.6 keV in Ref. [20j because of our choice 
of cosmological parameters. Besides, our sterile neutri- 
nos are hotter than those of Ref. J2(j , so the lower bound 



in our scenario should be even somewhat larger than 
2.9 keV. 

In the mass range in which sterile neutrinos can be 
part of the HDM, we can apply the bounds on the sum 
of the contributions of active and sterile neutrinos to the 
HDM density. The 3<r bound on the sum of the neutrino 
masses (see Fig. 2 of Ref. j^) is 

^m,+fm s <l.leV, i=l,2,3 (9) 

where rrii are the light neutrino masses. Combining Eq.[5] 
with an estimate of the light neutrino masses we obtain 
an upper limit on m s . If the neutrino mass spectrum 
is normal hierarchical, oscillation data impose the sum 
of the active neutrino masses to be about 0.05 eV. This 
provides the most conservative bound on m s : 

m s sin 2 29 < 0.1 d~ x (Tr/5 MeV) -3 eV, (10) 

which we plot for m s < 100 eV (vertically hatched region 
with thin lines in Figs. 1 and 2). 

The 3cr upper bound imposed by big bang nucleosyn- 
thesis (BBN) on any extra contribution to the energy 
density, parametrized as extra neutrino species, AiV„, is 
AN V < 0.73 (see Fig. 7 of Ref. |12), which translates into 
(horizontally hatched region with thick lines in Figs. 1 
and 2.) 

sin 2 29 < 5.6 d' 1 10~ 2 (Tr/5 MeV) -3 . (11) 

Let us turn now to experimental bounds. So far disap- 
pearance experiments have reported negative results |2.'l 
I24II25I]. Reactor-Pg experiments CHOOZ [23] and Bugcy 
|24j constrain the mixing angle relevant in v e — > D s con- 
version [2i|] to be sin 2 29 < 0.1, for m s > 1.7 eV. For 
smaller masses (down to m s ~ few 10 _1 eV) the bound 
is stronger by a factor of 2-5 |24j (darkest region in Fig. 
1). Accelerator-^ disappearance experiments |25| im- 
pose the mixing angle which controls — ► v s oscillations 
to be sin 2 29 < 0.02, for 13.8 eV < m s < 17.9 eV. Less 
stringent bounds apply for other values of m s (darkest 
region in Fi g. 2). F uture experiments looking for D e |2fij 
and v^, |27l l2Sj disappearance might explore part of 
the now allowed parameter space, e.g., ^-factories may 
reach sin 2 29 ~ 10~ 3 . 

Appearance experiments searching for v a — > vp oscil- 
lations (a, = e, /i, t), are sensitive to the product of the 
mixing angles between v ai up and the mostly-sterile mass 
eigenstate. These experiments have reported no positive 
signal, except for the LSND experiment (2j| , which found 
evidence of — > v e conversion. MiniBooNE [27J will test 
this result. Let us notice that in our model, the ranges 
of m s and sin 2 29 required to explain the LSND data (in 
terms of neutrino oscillations), are cosmologically and as- 
trophysically allowed. The analysis of the data requires 
the mixing of, at least, four neutrinos, and therefore can- 
not be used to set bounds on the mixing angle sin 2 29 we 
have used to parametrize 2—v oscillations 10] . 

In the mass range of interest, /3-decay experiments 
searching for kinks in the energy spectra of the emitted 
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electron constrain the mixing angle between v e and the 
mostly-sterile neutrino mass eigenstate. Different nuclei 
have been used and negative results have been found so 
far (for a complete review see Ref. |3(J)- The limits are 
strongly mass dependent (darkest region in Fig. 1). 

If neutrinos are Majorana particles, the neutrinoless 
double beta (((3f3)o v ) -decay is allowed. The half-life time 
depends on the effective Majorana mass (m) (see e.g. 
Ref. [IJ). The contribution of the mostly-sterile neutrino 
is of the form (m) s = m s sin 6e % P" where sin 2 29 is the 
mixing parameter in v e — > v s conversions and (3 S is a Ma- 
jorana CP-violating phase. At present, the most strin- 
gent bound on |(m)| is |(m}| < 0.35 — 1.05 eV j3^], which 
conservatively translates into (dashed line in Fig. 1): 

m s sin 2 29 < 4 eV . (12) 

Possible cancellations in | (m) | among contributions due 
to different mass eigenstates would weaken this bound. 
Barring this possibility, the present and future {P/3)q^- 
decay experiments (see, e.g., Ref. |33jp will probe part 
of the cosmologically relevant region, possibly up to 
m s sin 2 29 ~ 0.05 eV, in which v s could be an impor- 
tant part of the dark matter or produce pulsar kicks. 
The sterile neutrinos with sin 2 29 « 0.1 - 0.01 would 



have the cross section required in Ref. [34J to separate 
atmospheric showers produced by these neutrinos from 
showers generated by active neutrinos, in future experi- 
ments such as EUSO and OWL (35|, by using the Earth 
as a filter. The required flux of ultra-high energy neutri- 
nos would be very large. 

We presented here a scenario, based on a low reheating 
temperature at the end of inflation, in which the coupling 
of sterile neutrinos to active neutrinos can be as large as 
experimental bounds permit. For example, the sterile 
neutrino required by the LSND result does not have any 
cosmological or astrophysical problem. 

The experimental discovery of a sterile neutrino in 
the region of m s — sin 2 29 opened up in this paper, 
would require an unusual cosmology, such as one with 
a low reheating temperature as presented here. In this 
case, baryon asymmetry might be generated through the 
Affleck-Dine mechanism |36|, and the bulk of the dark 
matter (if not made of sterile neutrinos) should consist 
of other non-thermally produced particles. 
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